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ABSTRACT 

In the use and interpretation of log TV-log S distributions for gamma-ray 
bursts, burst peak flux has typically been used for S. We consider here the use 
of the ffuence as a measure of S, which may be a more appropriate quantity 
than the peak flux in such highly variable sources. We demonstrate how using 
the BATSE trigger data we can determine the selection effects on fluence. 
Then using techniques developed elsewhere to account for the important 
threshold effects and correlations. Applying the appropriate corrections to 
the distributions, we obtain a fluence distribution which shows a somewhat 
sharper break than the peak flux distribution, implying a possibly narrower 
fluence luminosity distribution. If bursts are at cosmological distances, these 
observations together indicate that evolution of the luminosity function is 
required. 

Subject headings: gamma rays: bursts 
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1. Introduction 

Because of their transient nature and lack of counterparts, the distance to Gamma-Ray 
Bursts (GRBs) is not known. Statistical distributions have therefore been used as indirect 
means of estimating their distances. The original discovery by the Burst and Transient 
Source Experiment (BATSE) team (Meegan et al. 1992), and the subsequent confirmations 
(Meegan 1996) that GRBs are isotropically distributed in the sky has gradually strengthened 
cosmological interpretations of GRB sources. If bursts are indeed of a cosmological origin, 
then the so-called logiV-logS' distribution can be used to constrain the model parameters. 

Under the assumption that bursts are distributed homogeneously and isotropically in 
a static Euclidean space (HISE), the logarithmic slope of the logiV-logS' distribution is 
expected to be -3/2. In cosmological scenarios, the observed deviations from this slope are 
due to the breakdown of the last two conditions in HISE. The degree of deviation depends 
on not only the geometry and expansion rate of space but also the shape and evolution of 
the luminosity function. This latter effect is less significant for narrow and slowly evolving 
luminosity functions, but if the range of luminosities is larger than the observed range 
of S, it will obscure the cosmological effects. Therefore the choice of the parameter S is 
important in the interpretation of the log TV-log S results. 

For steady sources, this choice is obvious; S would be represented by the steady and 
well-defined photon or energy flux. However, for highly transient sources such as GRBs, 
the measured flux depends on the time scale over which the burst is observed. Ideally, one 
would like the flux measure to be instantaneous, defined such that the observational time 
scale for the accumulation of photons is smaller than the intrinsic variation time scale of 
the source. In practice, observational measures such as the "peak flux" have been averaged 
over a time interval At which might not be small compared to the intrinsic time scale. Use 
of this version of the "peak flux" for the parameter S in the log iV-log S distributions can 
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lead to ambiguous interpretations, and a number of complicated steps have to be taken in 
order to extract the instantaneous peak flux distribution from the data (see Lee & Petrosian 
1996a, hereafter LP). 

These complications can be avoided by using a different observational measure for 
S. One such measure is the fluence J 7 , which we define to be the total (time-integrated) 
radiant energy per unit area within the BATSE trigger range of 50-300 keV. The fluence 
distribution may prove to be a more useful tool in cosmological studies than the peak 
flux distribution if the time-integrated luminosity has a narrower intrinsic dispersion and 
undergoes less evolution than the peak luminosity In fact, given the wide dispersion in 
the durations of GRBs, it is difficult to justify why the intensity of the highest spike in a 
bursting source should have a narrow distribution. Perhaps the total energy released has 
a narrower distribution than that of the peak luminosity. For example, it is more likely 
that the total energy released in a compact object merger might be a more appropriate 
"standard candle" than the rate of energy release in such a merger. Furthermore, in cosmic 
fireball models which presumably describe the physics of these mergers, the relationship 
between the luminosity observed in the detector's rest frame to that emitted by the source 
depends strongly on the bulk Lorentz factor of the expanding shell (Meszaros Sz Rees 1993, 
1994; Madras & Fenimore 1996), which is not likely to be universal for all bursts. 

In this paper we investigate the log iV-log T relation for BATSE GRBs. In the next 
section we describe how a bias-free fluence distribution can be obtained from instruments 
such as BATSE which trigger on some average flux value. In §H we discuss our results. 



2. Analysis 
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2.1. Obtaining the Fluence Limit 

We use the publicly available BATSE 3B catalog data, which provides the maximum 
and minimum photon counts C max and C min for the three trigger time intervals 
At = 64, 256, 1024 ms. The most fundamental selection effect is that C max must exceed the 
threshold C m i n . We are interested in the more physically meaningful fluences and fluxes. 
The BATSE catalog also gives values of the total energy fluence T within the 50-300 keV 
range, along with three measures of the average peak flux. The average peak flux is given by 
fp — C max /(A e ff(9,4>)At), where A e ff is the effective detector area (including the spectral 
response) of the instrument in the direction (6,4>). 

Given fp or JF, we ask what would have been the threshold for detection of a burst 
for any of these quantities. It is easy to see that a burst with average peak flux fp 
coming from a direction 6 and <fi would trigger BATSE if fp > fu m = C m i n / (A e ff(9, <p)At). 
Otherwise C max would be less than C min . Therefore, the threshold on the average peak 
flux is fn m = fpC min /C max . The same relation also exists for the fluence and its limit, as 
can be seen schematically by examining Figure |I[ A burst with observed fluence JF (or fp) 
and a particular pulse profile, spectrum, etc. would have been undetected if its fluence (or 
flux) was lowered by a factor of C max /C m i n , because then its peak counts C max would have 
been less than the limiting counts C m i n . Clearly then Tu m = FC m i n /C max , as long as the 
background count rate does not vary significantly throughout the duration of the burst. In 
summary, 

C m ax fp -F 

Cmin flim -F Mm 

The last equality is approximate because of the implicit conversion between photon counts 
and energy. This will be a good approximation if the source spectra do not change 
drastically throughout the duration of the burst. 

Note that all of the details of the pulse profile, spectrum, and instrumental response are 
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hidden in the ratio. The problem of obtaining the distribution of any of the quantities in 
equation ([]]) can be described generally as the problem of obtaining the distribution of some 
variable x subject to the condition that Xi > x^um for each data point i. A general solution 
to this problem was described by Petrosian (1993) and has been extensively discussed in 
LP. It can be seen that extracting the ffuence distribution is in principle no different from 
extracting the peak flux distribution. It should also be noted that the conclusions drawn 
from the distribution of V/V max or its average are unchanged no matter which of these 
properties is used as a measure of distance. 



2.2. Fluence Limit Interpretation 

Although the extraction of the fluence distribution is computationally straightforward, 
the interpretation of Tu m differs from that of C m i n and hence deserves some explanation. 
C m i n is simply a threshold that depends only on the background count rate and is a variable 
independent of the physical burst properties. For bursts which have durations T <C At it 
is clear that T = C max (his) /A e ff(9,4>) and Tu m = C m i n (his)/A e ff(0,(f>) are independent, 
where (hv) is the average photon energy in the 50-300 keV range. For long duration bursts 
with T > At, we have approximately T oc fpT(hu), so that from equation ([!]) we obtain 

T lim = (hu)T ° mi \ , (2) 

indicating that if as expected {hi;), A e ff(9,4>), and C min are essentially independent of J 7 , 
the fluence limit is approximately proportional to the duration. 



2.3. Fluence Distributions 

We now examine all bursts for which a fluence measurement and a value of C max jC^ 
exists. The values of C max /C m i n are known for three time scales: 64 ms; 256 ms; and 
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1024 ms. We use the 1024 ms values of C max /C min because this time scale is the most 
sensitive and allows us to use the largest number of bursts, but we also note that the other 
time scales give essentially identical results. Note that unlike in the determination of the 
distribution of instantaneous peak flux (see LP), where a correction for the short duration 
bias based on some observational duration measure was necessary, there is no need to have 
a separate measure of duration to determine the fluence limit. This increases the number of 
bursts available from 514 to 555 for the 1024 ms trigger. It may be argued that bursts with 
no well-defined durations may have data gaps or some other problems which would make 
the fluence measurements unreliable. As it turns out, the resulting fluence distributions are 
insensitive to whether or not we include the extra 41 bursts. 

The bivariate distribution of JF and Tu m is shown in the top panel of Figure |2|. 
Obviously because of the data truncation due to the variation in J-u m , simply binning 
the fluences to get a distribution will result in a biased distribution. As explained by 
Petrosian (1993), a nonparametric method exists to obtain a single variable distribution 
from a truncated bivariate distribution. This method amounts to using information from 
untruncated regions to estimate the data that was missed due to the truncation. Clearly, 
this is only possible provided the variables are uncorrelated. The first thing we must do is 
test the data for a correlation between T and Tu m . Using the correlation test designed for 
use on truncated data (Efron & Petrosian 1992), we find that the probability that the data 
are uncorrelated is 2.3 x 10~ 5 . As discussed in § [2.2| , since Tu m is approximately proportional 
to the duration, a correlation test involving T and Tu m effectively tests the correlation 
between fluence and duration. The results indicate that the fluence and the duration are 
positively correlated with each other. If the fluence is a good measure of distance, this 
result seems to be in the opposite sense of that expected from cosmological time dilation 
(c.f. Norris 1995). However, there are a number of factors which could complicate this 
interpretation (see Lee & Petrosian 1996b). 
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In order to go further, one must resort to parameterizing the correlation. As we have 
done before (Lee, Petrosian, & McTiernan 1993; 1995; LP), we use a simple power law 
parameterization. Briefly, we transform Tu m into T[ im = J-'umJ 7 ' and vary a until the 
correlation between T and T[ im disappears. This requirement gives a well-defined value for 
a. We find a = 0.22 ± 0.07, with the error interval indicating the ±90% confidence limits 
on a. The data truncation boundaries are transformed accordingly. The resulting bivariate 
distribution, which now contains uncorrelated variables, is shown in the bottom panel of 
Figure ^ which can then be readily integrated over T[ im with the methods described in LP. 
Using this technique, we obtain the cumulative and differential distributions of J 7 , along 
with the logarithmic slope of the cumulative distribution as a function of T (Fig. |3|). The 
qualitative shapes of these distributions persist even if we use different samples of bursts 
corresponding to higher fluence limits. We justify the particular parameterization chosen 
here (power law) by noting that similar values of a are found for data subsets chosen from 
various ranges in J 7 , and in any case the dispersion in the data is such that there would be 
little justification for a more complicated parameterization. 

Dividing our best estimate of the differential distribution n{T) by what would have 
been obtained without consideration of the truncation and correlation gives the trigger 
efficiency as a function of fluence, which is plotted in Figure [| along with the ratio of 
the observed number of bursts N obs (> J 7 ) to the total number of bursts iV(> J 7 ) greater 
than a given fluence J 7 . Our derived efficiency can be compared to the results of in't 
Zand & Fenimore (1996) and Bloom, Fenimore, & in't Zand (1996), who utilize a quite 
different approach. Rather than starting with the data and working backwards through 
the selection effects to derive the distributions, in't Zand & Fenimore use Monte Carlo 
simulations of a sample of bursts with a distribution of temporal and spectral shapes as 
observed by BATSE. Then assuming a cosmological origin for the sources, they predict the 
trigger efficiency of BATSE as a function of fluence. In contrast, our "efficiency" is purely 
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empirical and involves no model assumptions. Using these trigger efficiencies, Bloom et al. 
(1996) correct the observed fruence distributions to obtain the true distributions. Their 
results are qualitatively similar to ours, although their log iV-log .T 7 curve has a very steep 
upturn at low fluences which is not evidenced in our curves. This difference may arise as 
a result of their assumption that any correlations between burst characteristics are solely 
a result of cosmological effects such as time dilation and redshifting. However, it would be 
very difficult to reconcile the strong positive correlation that we found between fluence and 
duration with any non-evolving cosmological population of bursts. 



3. Discussion and Conclusions 

We have described a robust method of accounting for the selection biases on the 
detection of BATSE GRBs based on their fluences, independent of duration or spectral 
measures. Using methods described in our previous publications we have obtained the 
variation of the cumulative distribution, differential distribution, and logarithmic slope 
of the distribution as a function of fluence. The results shown in Figure |3] reveal that 
the fluence distribution appears to show a sharp break from slope -3/2 to about -1/2 
at T ~ 10~ 5 erg cm -2 . The analogous peak flux distributions (see Fig. 7 of LP for an 
example) show a slightly more gradual transition in slope. Our interpretation of this result 
would be that the time-integrated luminosity has a narrower distribution than the peak 
luminosity. Therefore the fluence may be a better indicator of the burst distance than the 
peak flux. The lack of consistency with the time dilation effect would imply that either the 
burst sources are not cosmological or models more complicated than the simple no-evolution 
model are necessary. 

Ignoring the inconsistency for the moment, fits of the fluence distribution to very 
simple cosmological models (constant comoving density, no luminosity evolution, energy 
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spectra either power laws or as in Band et al. (1993), density parameter ft — or 1, 
H = 75 km s _1 Mpc -1 ) give sources of total radiant energy £ ~ 10 52 ergs or £ ~ 10 51 ergs 
for Q = and Q = 1 models, respectively. Note that uncertainties in the spectral form can 
be absorbed into uncertainties in the luminosity evolution, and in any case the data are not 
sensitive enough to definitively distinguish among the models. A difference between these 
models and those involving peak fluxes is that the inferred maximum redshifts are greater 
{zmax ~ 3 for the models involving power law spectra and z max ps 5 for the models utilizing 
the Band spectral form), a result also noted by Bloom et al. (1996). However, adding in the 
evolution necessary for agreement with the time dilation results would reduce these inferred 
maximum redshifts. 

We acknowledge W. Azzam and G. Pendleton for useful discussions, and we thank the 
anonymous referee for comments and suggestions which led to an improved paper. This 
work was funded by NASA grants NAGW 2290 and NAG-5 2733. 
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Fig. 1. — Schematic diagram of burst selection. The larger histogram represents the light 
curve of a typical burst. Neglecting spectral variation effects, the fluence is proportional 
to the time integral of the curve. The shaded bin represents the peak flux of the burst 
(integrated over the trigger time), and the dotted line represents the average value of the 
limiting flux, below which the burst would not have triggered. The smaller histogram shows 
the light curve scaled down by a factor of fp/ fu m (or equivalently C max /C m i n ). No matter 
what the light curve looks like, it can be seen that if C max /C m i n is constant throughout 
the burst, then the limiting fluence for the burst is simply given by the fluence scaled by 

C m ax/ C m i n . 
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Fig. 2. — top panel: Bivariate distribution of JF and Tu m . The diagonal line indicates the 
selection criterion T > Tu m . bottom panel: Bivariate distribution of T and T' Um = Tu m T~ a : , 
with a = 0.23. The diagonal line indicates the selection criterion T > T' Um . 
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Fig. 3. — top panel: Cumulative distribution N^J 7 )^ 1 - 5 of fluence, showing the deviation 
from HISE. The solid histogram is the distribution that would have been obtained without 
consideration of selection effects or correlations. The dotted histogram is the distribution 
that would have been obtained by accounting for selection effects but neglecting correlation. 
The solid histogram is the distribution obtained when accounting for both effects. The 
dot-dashed line indicates the HISE prediction of logarithmic slope -1.5. middle panel: 
The differential distribution multiplied by JF 5 / 2 , which shows the deviation from the HISE 
prediction (dot-dashed line), bottom panel: Logarithmic slope of N as a function of T . The 
heavy solid lines show the 90% confidence limits on the transformation used to remove the 
correlation. 
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Fig. 4. — The lower two histograms show the trigger efficiency rio^F) /n^F) as a function 
of fluence, while the upper two histograms show the ratio of the cumulative number of bursts 
observed N obs (> T) to the true cumulative number of bursts iV(> J 7 ). The dashed lines 
depict our best estimates of the trigger efficiency and ratio of observed to total number of 
bursts, accounting for both data truncation and correlation. The dotted lines depict the 
efficiency or ratio without accounting for the effects of correlation. 



